OF1.7中的p |
您所在的位置:网站首页 › p等于肉gh h的单位 › OF1.7中的p |
翻译自:CFD-online 帖子地址:http://www.cfd-online.com/Forums/openfoam-solving/80454-p_rgh-1-7-a.html stawrogin: 亲爱的Foamer们 我对于OpenFOAM当中新的压力p_rgh有点困惑。它是怎样定义的?我有点担心的设置好的边界条件。比如如果我想要在出口就设置一个常数的"压力",这应该如何定义?如果定义p_rgh = 0,是否意味着在这个算例中pressure + (rho)*g*z是一个固定值? 我也很惊奇为什么在一些求解器中同时存在p—文件和p_rgh文件。P—文件仅仅用在后处理吗或者它也有可能在这里定义了边界条件? 先谢过 Stawrogin juho: 在1.7.x的buoyantPimpleFoam求解器中,p和p_rgh文件都被读取。P被用在热力学模型中。 在createFields中,p_rgh被定义为: p_rgh = p - rho*gh; 因此这个压力没有静压并且被从p文件中的压力场开始初始化。 在pEqn.H头文件中,压力方程被p_rgh来写和求解,因此p_rgh边界条件对于压力求解很重要。在解算完压力后,p通过下式计算: p = p_rgh + rho*gh; stawrogin: 亲爱的Juho 非常感谢你的帮助,你的回答对于我来说很清楚了。 Stawrogin deniggo: 你好,Juho p_rgh似乎是动压。我使用pimpleFoam求解器和pisoFoam求解器。是否可以调整这些代码获得p_rgh? 谢谢 Nico Angela Wang: 你好,Nico p_rgh代表动压,是否意味着p_rgh=1/2*rho*U^2? 我检查了案例中damBreak的算例,为什么p_rgh的值比p的值大?这是否意味着rho*g*h是负数? 我非常困惑。 Quote: Originally
Posted by deniggo Hi Juho, p_rgh seems to be the dynamic pressure. I work with pimpleFoam and pisoFoam. Is it possible to adjust these codes to obtain p_rgh? Thanks, Nico kjmaki: 你好,Angela p_rgh是p-rho*g*h或动压。它在1.5版本中被叫做pd,它在版本1.6中是为了求解总压,在版本1.7中仅仅是为了求解动压。 Kevin deniggo: 你好,Kevin 但是为什么p_rgh比p高?一个负的静压(rho*g*h)怎么可能? 不应该是1/2*rho*U^2(动压定义)导出p_rgh?贴出计算后U—场的分布。 nico nbadano: 亲爱的朋友们, 够古怪了,gh在OpenFOAM中被定义: gh = g & mesh.C() 换句话说,gh是g向量和单元中心位置向量的点积。g常常被定义为(0 0 -9.81),gh的值常常正z是负值! 结果的导致: p_rgh = p-rho * gh 如果z>0,p_rgh比p大 希望这个有帮助 祝好 Nicolas Angela Wang: 我的问题解决了,p_rgh是动压= p- rho*g*h Logan Page: 关于这个一个问题: 使用buoyantBoussinesqSimpleFoam求解器,p和p_rgh是(m^2)/(s^2)[比如流动压力] 很明显的是这样定义太简单了:(p_rgh/rho)=(p/rho)-g*h 我的问题是当通过划分,密度如何使用。是rhok = 1 - beta (T - Tref)?? 谢谢 kerim: 朋友们好 我将要继续Nico的想法。对于我来说,p_rgh不是动压。我将层流溃坝算例的结束时间改为500。那时初始流后将固定。因此在t=500时没有动压。但是interfoam求解器没有给定一个非0的p_rgh压力。这就是为什么我认为p_rgh不是动压。谁能告诉我发生了什么吗?或者我哪里错了吗? Kerim nbadano 你好,Kerim p_rgh不是真正的动压,特别是在整个计算域内密度变化的多相流。是压力场和rho*g*z场。我认为这也是这个场在1.5版本没有被叫为pd的原因。 为了增加困惑,记住rho*g*z也不是真的静压!!甚至密度是常数(在z=0和大气压力p=0平面之间乘以rho*g)。 这里有一个p_rgh, rgh和p作为静压条件的一个草图(完全没有移动)。希望有帮助。 Cristiano: 这里的朋友你们好 我正在模拟一个在大气条件下测试过的燃烧器并且我有点困惑。 我应该如何为p和p_rgh设置合适的边界条件? Inlet --> Outlet-> Wall---> 实在谢谢你 Cristiano linch: Quote: @nbadano I agree to your sketch: if the p-field is smooth, the h-field is also smooth, the g-field is uniform and there is a jump in the rho-field at h=h2, so the must be a jump in thep_rgh field. But if you set up your example in interFoam you'll notice, that the p_rghfield is uniform also. Though the p_rgh definition might be a littlte different. Well, theoretically you'll get an uniform field, if you divide p_rgh by rho, but in this case you wouldn't get the p-dimension [kg*m^-1*s^-2]. So somehow the hydrostatic pressure is being substracted and not the simple product of rho*g*z. 对不起,我错了。P_rgh不是均匀的。 christianfrias: 尽管这个建议很简单,但还是值得说。 你有一个成功运行OpenFOAM1.6版本的算例并且想要运行在OpenFOAM1.7平台。记住你仅仅需要重命名你的p文件为p_rgh(在你的transportProperties 和fvSchemes文件里面改变一些其他参数)。然后将g的值从(0 0 -9.81)改为(0 0 0)。这个算例应该可以和OpenFOAM1.6一样运行,如果他们能顺利运行在OpenFOAM1.6下,你也无需担心你的边界条件是否适用于OpenFOAM1.7。 祝好 Christian F 为了理解不同压力,请看伯努利: 动压 --> 1/2*rho*v^2 液压--> rho*g*h 静压 --> p 1/2*rho*v^2 + rho*g*h + p = Constant 从OpenFOAM的站点,p_rgh=p-rho*g*h。 因此,p_rgh是静压减去液压,在任意高度 我不知道静压是从这个讨论的什么地方出来的。我认为p_rgh没有任何意思。动压是运动流体的压力并且在流体速度为0时,它将和静压转换。能量守恒回到伯努利。 很高兴知道: 为了让我的VOF模型运行,我将fvSolution文件放置在PISO求解器下。 pRefPoint (0.0 0.0 0.0); pRefValue 1e5; 成功! deniggo: 大家好 谢谢King清楚的解释。 为了解开压力疑惑:后处理工具"ptot"计算了每个时间步的总压(静压+动压): http://www.cfd-online.com/Forums/openfoam/74978-where-define-fields-write.html 祝好 Nico fasfcastro: 我认为p_rgh是扰动压强(从一个流体力学平衡状态)实际被用于模拟动量方程。 Oke'e: 为什么在bouyantBoussinesqPimpleFoam求解器和bouyantBoussinesqSimpleFoam求解器算例中p_rgh的单位是m^2s^-2? 但是在(多相流)interFoam算例中的单位是kgm^-1s^-2。 Andrea_85: 因为在前两个算例中除以了密度。为了获得"真实"压力必须乘以密度[kg/m^3]。 祝好 andrea MOHAMMAD67: 亲爱的朋友们 你们好 这里已经有了许多关于p和p_rgh的讨论。但是自己没有从这些讨论中得出最后的结果。我想知道在多相流中为了和实验对比,我应该选择什么压力?如你所知,我们在实验室里只能得到液压。 p_rhg的值是否依赖于初始的计算域? christianfrias: 我看既然p_rgh不是动压也就意味着p不是全压(如我们所想)实际它是静压。为了计算全压,你可以使用ptot(它将会计算1/2*rho*U^2 [动压] + p [静压)])。因此,为了将OpenFOAM的结果与实验数据对比,我将使用OpenFOAM中代表静压的p文件。 sharonyue: Quote: Originally
Posted by christianfrias I see that since p_rgh is not the dynamic pressure it means that p is not the total pressure (as we could think) and is actually the static pressure. To calculate the total pressure you can use ptot (this will calculate 1/2*rho*U^2 [dynamic pressure] + p [static pressure)]). So, to compare experimental data with the results from OpenFOAM I will use the p file which is the static pressure in OpenFOAM. 因此p(静压)加上1/2*rho*U^2(动压)常常在实验中被测量?为什么我记得在一个方盒或者其他什么东西,静压被一个压力表测量。 vonboett: Quote: Originally
Posted by sharonyue so is the p(static pressure) plus 1/2*rho*U^2(dynamic pressure) usually measured in experiment? Why I remember usually in a tube or something ,static pressure is measured in a pressure gauge 依据你测量的什么。水压力水平传感器(和压强计)眼下等于液压关于传感器,因此水平液压传感器并且你需要和OF中的p进行对比。如果你测量的是流动方向阻力的增加等。这个就可能很难直接对比测量和模拟了。 baedmaen: 大家好 我现在对于多相流的压力很困惑。有人能帮我把定性压力特征和附件中的文件关联起来吗? 第一幅图: 静止的两相系统(水和空气);我认为nbadano已经在2010年12月3日通过帖子给出了答案。 第二幅图: 气泡上升的表面,截图 第三附图: 一滴水滴入水的表面 非常感谢 附件: Dream: 我们设置了分离器模型并且模拟了多相流。在我们运行了setFields,下面的信息出现了: Setting field default values Setting internal values of volScalarField alphaair Setting internal values of volScalarField alphawater
--> FOAM FATAL IO ERROR: Cannot find patchField entry for wall
file: /home/dell/OpenFOAM/damBreak4phaseFinelmxlmx/0/alphawater.boundaryField from line 34 to line 56.
From function GeometricField::GeometricBoundaryField::readField(const DimensionedField&, const dictionary&) in file /home/opencfd/OpenFOAM/OpenFOAM-2.2.2/src/OpenFOAM/lnInclude/GeometricBoundaryField.C at line 206.
FOAM exiting -----------------------------------alphawater------------------------------------------- FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object alphawater; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 1;
boundaryField { Wall { type zeroGradient; //value uniform 0; } new_new_out { type fixedValue;//; value uniform 1; } new_out { type fixedValue; value uniform 1; } out { type fixedValue;//outletInlet; value uniform 1; } in { type outletOutlet; outletValue uniform 0; value uniform 0; } }
可能我们的边界条件设置不正确,但是我们不知道如何修改他们。 非常感谢帮助 Andrea_85: 你好 用 Wall { type zeroGradient; //value uniform 0; } 替换 wall { type zeroGradient; //value uniform 0; } 祝好 andrea vonboett: Quote: Originally
Posted by baedmaen Hi all, now i am confused about multiphase pressure. Could anybody put the corresponding qualitative pressure profiles into the file attached? First figure: stationary two phase system (water and air); i think nbadano already posted the answer at December 3, 2010, via this thread Second figure: bubble rising to surface; snap shot Third picture; first contact of a water drop with water surface Thank you very much! 好的,p_rgh依赖于你的坐标系统,边界条件和最后的pRefValue和pRefPoint。然而,如果你有一些像大气pressure=0Pa在其下为空气通过指定大气边界或者pRefValue = 0 Pa,你的水密度和g被正确指定,你会得到一个的很好的从0到表面压力p三角液压。我在纳维斯托克斯方程方程中使用p_rgh但是我使用p作为一个比较结果。我的附件是两张设置大气环境为0Pa,3D计算域从底部z=2m到顶部z=3m,一个气泡在z=2.25m处,同时一个半径为0.01m的水滴在z=2.51m处。自由面在2.5m处。这个截图是在开始模拟T=0.01s抓取的。在各个方向上的单元尺寸为5mm。注意p_rgh的密度影响,20468Pa=2.51m*g*rho_water。 Dream: Quote: Originally
Posted by Andrea_85 Hi, replace Wall { type zeroGradient; //value uniform 0; } with wall { type zeroGradient; //value uniform 0; } Best, andrea 非常感谢! 我们解决了这个问题,原因在于我们大写了单词壁面的首字母W John Handel Kennedy: 大家好, 我打算模拟直管中的换热 我使用buoyantBoussinesqSimpleFoam求解器 我设置g和beta为0 我的壁面温度是373K并且入口温度为293K 入口的速度为1m/s 管的直径为1m并且为了让雷诺数达到100设置nu的值为0.01 层流普朗特数为1.5 我用simpleFoam求解器得到了充分发展的流动。例如速度增加到2m/s 然而我用buoyantBoussinesqSimpleFoam求解器不能得到相同的速度分布,速度朝出口降低。 我应该如何解决这个问题? 我应该在alpha_t和p_rgh文件中指定什么吗? 祝好 John |
CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3 |